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Abstract. We theoretically propose an alternative method to realize a deterministic 
SWAP gate using shortcuts to adiabatic passage based on the approach of Lewis- 
Riesenfeld invariants in cavity quantum electronic dynamics (QED). By combining 
Lewis-Riesenfeld invariants with quantum Zeno dynamics, the SWAP gate can be 
achieved deterministically. The strict numerical results show that our scheme is a fast 
and robust approach to achieve SWAP gate. 
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1. Introduction 

With the development of quantum science, the quantum computer replaces the universal 
computer just around the corner. The quantum logic gate is an indispensable part of 
the quantum computation. In quantum computing especially the calculation model of 
quantum circuit, the quantum logic gate is a basic operation of qubits quantum circuit. 
Recently, a number of schemes have been proposed to perform quantum logic gates 
and quantum information processing using optical devices [1], quantum dot [2, [3], QED 
systems [4], ion trap and superconducting devices la m Ei m i- As we know that 
a universal set of quantum operations can be constructed by a series of single- and 
two-qubit gates. Even though a multi-qubits gate could be decomposed into one- and 
two-qubit gates, it will be extremely complex and time-consuming in practical operation, 
that the arising decoherence will destroy the quantum system eventually. An efficient 
method of quantum information processing requires the whole operation is robust 
against the decoherence, easily prepared and measured. Adiabatic passage technique 
provides the robustness of the method against small variations of field parameters, and 
the decoherence caused by spontaneous emission can be avoided if the dynamics follows 
dark states, i.e. states without components on lossy excited states. However, that is 
not precise. If the controlling parameters do not change slowly enough, the transitions 
between different time-dependent instantaneous eigenstates may still exist, that the 
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fidelity of the evolved state with respect to the target one must be reduced. In other 
words, the adiabatic passage technique usually needs a long process ms. If the required 
evolution time is too long, the speed of the system evolution will be slowed down, that 
the dissipation caused by decoherence, noise, and losses would destroy the expected 
dynamics finally. 

Shortcuts to adiabatic passage is a promising technique for quantum information 
processing which actually fight against the decoherence, noise, or losses that are 
accumulated during a long operation time. Thus, a variety of schemes have been 
proposed to construct shortcuts to adiabatic passage in both theories and experiment 

The shortcuts to adiabatic passage of logical gates 
operation have been presented, too. Chen et al [19] proposed a scheme of shortcuts to 
adiabatic passage for performing a tt phase gate. And we have proposed a scheme of 
shortcut to adiabatic passage for constructing the multiqubit controlled phase gate m 
In this paper, we effectively combine the advantages of shortcuts to adiabatic 
passage and quantum Zeno dynamic (QZD) [23] El] to implement a SWAP gate. It dose 
not need the composition of element gates, but directly implements the deterministic 
SWAP gate through designing resonant laser pulses by the invariant-based inverse 
engineering. The logical SWAP gates in our scheme can be performed in a much shorter 
time than that based on adiabatic passage technique, and it is very robustness to the 
decoherence caused by atomic spontaneous emission and cavity decay. 

This paper is structured as follows: In section 2, we give a brief description 
of the preliminary theory about Lewis-Riesenfcld invariants and QZD. In section 3, 
we effectively combine the shortcuts to adiabatic passage and QZD to implement 
the deterministic SWAP gate. Section 4 shows the numerical simulation results and 
feasibility analysis. The conclusion appears in section 5. 

2. Preliminary theory 

2.1. Lewis-Riesenfeld invariants 

We first make a brief introduction of the Lewis-Riesenfcld invariants theory g5] m\. 
Considering a system which is governed by a time-dependent Hamiltonian H(t), and 
we can seek the time-dependent Hermitian invariants I ( t ) that is related to the original 
Hamiltonian H{t) to satisfy 

ih^ = [H(t),I(t)]. (1) 

Obviously, its expectation values remain constant all the time, and drives the system 
state evolve along the initial eigenstate of /(t). For the time-dependent Schrodinger 
equation ih klALll — H(t)), the solution can be expressed by the superposition of 
dynamical modes |4> n (t)) of the invariants /(f) 

i m) = 


n 


( 2 ) 
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Figure 1 . (Color online) The schematic setup of SWAP gate implementation. The 
atoms are trapped in a single mode optical cavity, and each atom possesses five atomic 
levels. 


where n = 0,1,..., and C n is one of the time-independent amplitudes, a n is the Lewis- 
Riesenfcld phase. |<f> n (f)) is one of the orthonormal eigenvectors of the invariant I(t ) 
with the corresponding real eigenvalue A„, satisfying /(f)|<3> n (f)) = A n |<h n (t)). And the 
Lewis-Riesenfeld phase satisfies 

<*„(«)=t - H(t') i*„(t , )>. (3) 

2.2. Quantum Zeno dynamics 

The quantum Zeno effect was first proposed by Misra and Sudarshan, and it can be used 
to reduce the influence of the decoherence and dissipation by inhibiting the transition 
between quantum states through the frequent measurements EZf We consider a system 
which is governed by the Hamiltonian 

H k = H ohs + KH meas , (4) 

where H 0 bs is the Hamiltonian of the investigated quantum system and the H meas is the 
interaction Hamiltonian performing the measurement. K is a coupling constant, and 
when it satisfies K —> oo, the whole system is governed by the evolution operator 

U(t ) = exp [—it YjK\ n P n + P n H ohs P n )\, (5) 

n 

where P n is one of the eigenprojections of H meas with eigenvalues A n (-H meas = A n P n ). 

3. Shortcuts to adiabatic passage for the deterministic SWAP gate 

We consider a system in which the atoms are fixed inside an optical cavity as shown in 
figure 1. Each atom possesses three ground states |0), 11), |a) and two excited states |e), 
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| u). The transitions |0) B- \e) and |a) B- |e) are coupled to the laser pulses, with the 
corresponding Rabi frequencies f2 0 (f) and £l a (t), respectively. The transition |1) B- |e) 
is strongly conpled to the single mode cavity field with the coupling constant g. And 
the auxiliary excited state | u) is only used to implement the one-qubit operation. 

Now we show the scheme to realize a deterministic SWAP gate. The initial state 
|T 0 ) of the two atoms and the cavity field is defined as 

|*o> = a oo|00)As|0)c + «oi|01)ab|0)c 

+ Q!io|10)ab|0)c + q;ii|11)ab|0)c, (6) 

where \nm) ab\$)c denotes atoms A and B are in the state | nm) AB ( n,m = 0,1), and 
the cavity is in vacuum state |0)c- aw denotes the amplitude of the state \nm) AB \0)c, 
and satisfies the normalization condition. The SWAP gate makes the values of the two 
qubits exchange and then the output state becomes 

I'f) = aoo|00)As|0)c + o:oi|10)ab|0)c 

+ cqo|01)As|0)c + o:ii|11)ab|0)o (7) 

There are three steps for the implementation of SWAP gate. First, we transfer 
the population of |01 )ab|0)c to —|1 ci)ab|0)c completely with the help of laser pulses 
resonant with A atomic transition |0 )a B- \e) A and B atomic transition | o)b B- \e)s 
with the corresponding Rabi frequencies Oo a(£) and Q aB (t), respectively. After the 
interaction, the initial state becomes 

l^i) = «Oo|00)ab|0)c — c^oi 11®) as |0) c 

+ «io|10)ab|0)c + «h|11)ab|0)c- (8) 

Then the population of |10 )ab|0)c is completely transferred to — |o1)ab|0)c with 
the similar method with the help of laser pulses resonant with B atomic transition 
|0 )b -B- | e)s and A atomic transition |o)a -B- \e) A with the corresponding Rabi 
frequencies VL QB {t) and fl aA (t), respectively. And the state of the system becomes 

I HA) = «oo|00)a_b|0)c — otoi\la) A b\Q) c 

— atio\al) A b\Q) c + Q'ii|H)as|0)c'. (9) 

Finally the population of the single atom state | a) A (B) is transferred into — |0 )a(b) 
by the one-qubit operation with the help of the auxiliary excited state |w)a(b), and the 
corresponding laser pulses are resonant with A(B) atomic transition |0 )a(_b) -B |m)a(b) 
and \a) A (B) -B |m)a(B) with the corresponding Rabi frequencies ^' aA ^ B )(^) anc ^ A(B)(^)i 
respectively. As a result, the state of the system turns into 

1^3) = <aoo|00)ab|0)c + aoi|10)As|0)c 

+ «lo|01)al|0)c + ail|ll)As|0)c, (10) 

thus, the SWAP gate is obtained. 

In the following, we show the details of how to realize the SWAP gate by combining 
the shortcuts to adiabatic passage and QZD. For the first step, the input laser pulses are 
resonant with A atomic transition |0 )a •B \e) A and B atomic transition \cl)b •b |e)# with 
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the corresponding Rabi frequencies fl 0 a( 0 and Q aB (t), respectively, 
of the first step in the interaction picture can be written as (Ti = 1) 

The Hamiltonian 

H 1 = H l + H c , 

(ii) 

Hi = fioA(i)|e)A(0| + Cl aB (t)\e) B (a\ + H.c., 

(12) 

H c = g A a\e) A (l\ + g B a |e) B (l| + H.c., 

(13) 


here a is the annihilation operator, and 9a(b) is the coupling strength between cavity 
mode and the trapped atom. In the following, we will set Qa = 9b for simplicity. 
Since the lasers do not couple with the atomic state 1 1), that the states |11)ab|0)c and 
|10)ab|0)c ar e decoupled from the other states in this step, and will not participate the 
evolution any more. For the state |00)ab|0)c, the evolution subspace can be spanned 
by the basis vectors |v?i) = |00)ab|0)c, W 2 ) = |e0)Afi|0)c, and |<^ 3 ) = |10)ab|1)c- The 
Hamiltonian of this system reads 

H\ = fIoA(^)|e)A(0| + <7 a|c)a( 1| + H.c. (14) 

For the condition gA Hoa(^), the Hilbert space is split into three invariant subspaces 
r 1 = 1^1), r 2 = ^(-|^ 2 ) + 1^3)) and r 3 = ^-(|<p 2 ) + 1^3)) with the eigenvalues 
Ai = 0, A 2 = —gA, and A 3 = gA- The corresponding projections are P“ = |a)(a| 
(|a) G r n , n = 1,2,3). According to the QZD, U(t) = exp[-it^2 n (KX n P n + P n H ohs P n )], 
the effective Hamiltonian reduces to 

H eS = g A \e)A{l\ + Su|l)A(e|, (15) 

which indicates that the transition between |00)ab|0)c and |eO)AB|0)c is eliminated. For 
another state |01)ab|0)c, the system evolves in the subspace which is spanned by the 
basis vectors |0i) = |01)ab|0)c, \ fa) = |el)As|0)c, \ fa) = |11)ab|1)c, Ifa) = |le)AB|0)c, 
and |0 5 ) = |1o)ab|0)c- This evolution system is governed by the Hamiltonian 

H 2 = fI 0 A(t)|e)A(0| + 5u|e)A(l| 

+ Q aB (t)\e) B (a\ + g B \e) B (l\ + H.c. (16) 

By using the similar way to the above processes, the effective Hamiltonian of this system 
reads 

H'eft = ^=|/*)(^Oa(*)(0i| + ftai?(t)(05 I + H.c), (17) 

where \/i) = + \4>4))- To speed up the transition from l^) to — \(p 5 ) by the 

invariant-based inverse engineering, we must find out the invariant Hermitian operator 
I(t) which satisfies (H]). Since the effective Hamiltonian P[' cS possesses the SU(2) 
dynamical symmetry, the invariant I(t) can be given by j28j 

I(t) = —7=x(cosz/sin/3|/x)(0i| + cos v cos/3|/i)(05 
\/2 

+ i sini/|0 5 )(0i| + H.c), 


( 18 ) 
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where x is an arbitrary constant with units of frequency to keep /(f) with dimensions 
of energy. The time-dependent auxiliary parameters v and f3 satisfy the equations 


z> = -j=(Q 0A (t) cos f3 - Q aB (t)sin/d), 

(3 = —= tan u(fl aB {t) cos (3 + hloA^) sin [3). 
\/2 


(19) 


From (fl9l) the expressions of flou(f) and flos(f) can be derived as follow 
hloA(f) = V2($ cot v sin (3 + v cos (3 ), 

^ a B (f) = V^2(/3 cot v cos (3 — v sin (3) , (20) 

where the dot denotes a time derivative. The eigenstates of the invariant I (f) are 


|$o(*)> = cos i/ cos /310i ) — i sinz/l/i) — cos z/sin 0|0 5 ), 

|<l>-i-(t)) = —/= [(sin v cos /3 ± i sin/?)^!) + i cosz/|/i) 
v 2 

— (sin v sin/3 =|= i cos/3)|0 5 )], 


( 21 ) 


with the eugenvalues £ 0 = 0 and e± — ±1, respectively. The solution of the Schrodinger 
equation ih = //(f)|T(f)) can be written as the superposition of the eigenstates 

of /(f) 


n=0,± 


( 22 ) 


where a n (t) is the Lewis-Riesenfeld phase in (3), and C n is a time-independent 
amplitude. In order to get the final state — |0s), we choose the parameters appropriately 

7 rf 

<t) = c, m = W , (23) 


where e is a time-independent small value. From the precise calculation, we can easily 
obtain 


SMf) 

^ae(t) 


7r 

V^f/ 

7T 

v^f/ 


cot e sin 


cot e cos 


7Tf 

Tit 
2 f f 


(24) 


When t = tf, 


l^(f/)) — — sine sin a|0!) 

+ (—i sin e cos e + i sin e cos e cos ck) |/r) 

+ (—cos 2 e — sin 2 ecosa)|0 5 ), (25) 


where a = 7r/(2sine) = |a±|, and a± are the Lewis-Riesenfeld phases. As long as a 
satisfies the condition a = 2Nti(N = 1,2,3...), |T(fj)) = —10 5 ) = — |1o)ab|0)c can 
be achieved. Then the hrst step is realized successfully, and |Ti) = o;oo|00)ab|0)c — 
aoi|la)As|0)c + o?io|10)us|0)c + «h|11)as|0)c can be obtained. 

The second step is just similar to the hrst step. In this step, the laser pulses are 
resonant with B atomic transition |0)s <->■ \e) B and A atomic transition |a)u -B- |e)u 
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Figure 2. (Color online) The three steps for constructing the SWAP gate. For each 
step, the initial state is denoted by an empty circle and the final state is represented 
by a black circle. 


with the corresponding Rabi frequencies fl 0 s(t) = cot e sin and £l a A(t) = 

cot e cos respectively. In this step, the states |00 )ab|0)c, — |1o)ab|0)c and 
|11 )ab|0)c will be remained the same, while the population of the state |10 )ab|0)c 
is completely transferred to the state —|al)An|0)c, then |\P 2 ) = ctoo| 00 ) J 4 B| 0 )c — 
aoi|la)AB|0)c - aio|al)As|0)c + Q!ii|11)ab|0)c can be achieved. 

The third step is just a one-qubit operation with the help of the auxiliary state 
| u)a(b)- The input laser pulses are resonant with A(B) atomic transition |0 )a(.b) tA 
| u)a(b) and \a) A {B) tA | u)a(b) with the corresponding Rabi frequencies £l a A(B){t)' = 
cot e sin -jjy and &q A (b) (t) = yfy cot e cos respectively. The Hamiltonian in this 
step is given as follow 

H 3 = fi oiWI“)<(a| +^oiWI«)i(0| + H.c. (26) 

i=A,B 

For the initial atomic state |a)A(n), the population will finally transfer to the state 
— |0 )a(b) by using the similar way in the first step. Thus, we can obtain |T 3 ) = 
«oo|00)as|0)c + «oi|10)ab|0)c + aio|01) a i|0)c + O!ii|11)ab|0)c- That is the result of 
the SWAP gate, figure 2 represents these three steps for constructing the SWAP gate. 
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Qi (t) / go 





Figure 3. (Color online) (a) Temporal profile of the time dependence Rabi frequencies 
fli{t)/go versus t/go with f li(t) = f 2 0 a(£) (dash blue line), O a s(f) (solid blue line), 
flo s(t) (dash red line), f l a A(t) (solid red line), fl aA (solid green line), fl 0 A(B){t) 
(dash green line), (b) Time evolutions of the populations of the corresponding system 
states with the initial states 101). (c) Time evolutions of the populations of the 

corresponding system states with the initial state 110). The system parameters are 
set to be e = 0.25, gA = 9 b = 10go and tf = 20/g 0 . 



Figure 4. (Color online) The effect of atomic spontaneous emission 7 on the fidelity 
of the SWAP gate with different values of the cavity decay k. 
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4. Numerical simulations and feasibility analysis 

In the following, we make the numerical simulations to verify the validity of the SWAP 
gate, figure 3(a) shows the time-dependence laser pulse Q m (t)/g 0 as a function of t/g 0 
for a fixed value e = 0.25, gA = 9 b = 10go and tf = 20/go- With these parameters, the 
Zeno condition can be met well. The populations of the states |01)ab|0)c and |10)ab|0)c 
swap perfectly which can be seen from figure 3(a) and figure 3(b). Whether a scheme 
is available largely depends on the robustness to the loss and decoherence, so in the 
following, we consider the loss and decoherence effects on our SWAP gate operation. 
The corresponding master equation for the whole system density matrix p(t) has the 
following form: 


pit) = - i[H totah p(t)\ - ~[a ] ap{t) - 2 ap{t)a ] + p(t)a j a\ 



l=e,u k=0,l,a 


(27) 


where i/totai = Hi + H 2 + H 3 . n is the cavity decay rate, 7 a(b) is A(B) atomic 
spontaneous emission rate from the excited state = e,u\i = A, B) to the ground 
state | k)i(k = 0,1, a), respectively, ay*, = |/) (k |. For simplicity, we assume ja = 7 # = 7 
and the initial condition p( 0) = |To) (4 / o|. In figure 4, the fidelity of the SWAP 
gate is plotted versus the dimensionless parameter 7 jg with different values of n/g 0 
by numerically solving the master (1271) . From figure 4 we can see that the fidelity of our 
SWAP gate is higher than 97.6% even when the values of 7 and n are comparable to 
7 = O.lgA(B) = go an d k = gA{B) = 10go- It shows that the SWAP gate in our scheme 
is robust against decoherence due to cavity decay and atomic spontaneous emission. 

Now we give a brief analysis of the feasibility in experiment for this scheme. The 
scheme can be realized with trapped ions and nitrogen-vacancy color center in diamond 
[ 29] , cavity QED systems [30, 31) [32] or with impurity levels in a solid, such as Pr 3+ 
ions in YoSiOs crystal [33]. In recent experiments, the fabrication of various liigh- 
Q microcavities including whispering-gallery-mode cavities !3Tj. 35], micropost cavities, 
and one- or two-dimensional photonic-crystal microcavities [3.6, 37] also have been well 
designed. And the suitable parameters of toroidal microcavity system for strong¬ 
coupling cavity QED have been investigated and the parameters can be chosen as 
(g, hi, 7 )/ 27 t = (750,3.5,2.62) MHz [35]. In this case, 7 ~ 0.0034g,4( B ) = 0.034g 0 , 
k, ~ 0.0046gyi(s) = 0.046g 0 , and with these parameters the fidelity of SWAP gate 
can reach 99%. So our scheme is robust against both the cavity decay and atomic 
spontaneous emission and could be very promising and useful for quantum information 
processing. 
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5. Conclusion 

In summary, we have proposed a promising scheme to implement a SWAP gate through 
the shortcut to adiabatic passage and QZD instead of relying on the compositions of a 
large number of elementary gates. We also study the influences of atomic spontaneous 
emission and cavity decay on the fidelity through numerical simulation. The numerical 
simulation results denote that our scheme is very robust against the decoherence caused 
by atomic spontaneous emission and cavity decay, so it can be a more reliable choice in 
experiment. We believe that our scheme will be useful to realize quantum algorithms, 
such as Shor’s algorithm for prime factoring or Grover’s algorithm for database search. 
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